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PRELIMINARY SCALING LAWS FOR PLASMA CURRENT, ION KINETIC 
TEMPERATURE, AND PLASMA NUMBER DENSITY IN THE 
NASA LEWIS BUMPY TORUS PLASMA 
by J. Reece Roth 
Lewis Research Center 

ABSTRACT 

Parametric variation of independent variables which may affect the charac- 
teristics of the NASA Lewis Bumpy Torus plasma have identified those which 
have a significant effect on the plasma current, ion kinetic temperature, and 
plasma number density, and those which do not. Empirical power- law correla- 
tions of the plasma current, and the ion kinetic temperature and number density 
were obtained as functions of the potential applied to the midplane electrode rings, 
the background neutral gas pressure, and the magnetic field strength. Additional 
parameters studied include the type of gas, the polarity of the midplane electrode 
rings (and hence the direction of the radial electric field), the mode of plasma od- 
eration, and the method of measuring the plasma number density. No significant 
departures from the scaling laws appear to occur at the highest ion kinetic tem- 
peratures or number densities obtained to date. 

INTRODUCTION 

The characteristics and performance of the superconducting NASA Lewis 
Bumpy Torus Magnet Facility have been described elsewhere (refs, 1 and 2). 

The properties and behavior of the plasma created in this facility have been dis- 
cussed in a series of reports (refs, 3 to 12), The NASA Lewis Bumpy Torus ex- 
periment is characterized by three factors: 

(1) The plasma, magnetic field, and ion heating mechanism are operated in 
the steady state.! 

(2) The ion kinetic temperature is typically more than a. factor of 10 higher 
than the electron temperature, 

(3) The plasma is acted on by a combination of strong dc magnetic and elec- 
tric fields. 

The presence of strong radial electric fields acting on the plasma, which 
usually exceed values of 1 kilovolt per centimeter, not only is responsible for 
heating the ions to kinetic temperatures on the order of kilovolts, but also affect 
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the dynamics and radial transport of the plasma. Previous investigations (refs. 

13 to 16) have shown that, in common with Penning discharges and magnetron- like 
devices, the plasma forms rotating spokes which gyrate around the minor circum- 
ference of the plasma with velocities comparable to the E/B drift velocity; and 
the ions in these rotating spokes then form an energy reservoir which is thermal- 
ized to the high kinetic temperatures observed. The strong radial electric fields 
also have a profound effect on the radial transport of plasma, 

Prer -'ous investigations (refs. 3 and 10) have shown that the average particle 
resided ’me of the ions is virtually independent of the magnetic field over a 
factor of 10 variation in the magnetic field strength. This previous work indi- 
cated that the ion heating and radial transport processes in the NASA Lewis 
Bumpy Torus differ substantially from those in stellarators, tokamaks, and other 
conventional approaches which rely on pure magnetic containment of the plasma. 

It was therefore deemed of interest to investigate the parametric scaling laws for 
the plasma current and the ion kinetic temperature and number density in the 
Bumpy Torus plasma to see whether scaling laws could be obtained over a range 
of an order of magnitude or more in the independent variables which have been 
investigated thus far. Such empirically derived scaling laws might afford a clue 
to the physical processes responsible for the radial transport of plasma in this 
device, and also might confirm the ion heating mechanism which has been investi- 
gated both theoretically and experimentally on a prototype mirror machine facility 
(refs. 14 and 15), The data on which the correlations of this report are based 
were taken over a period of 2-| years, and amounted to nearly 800 individual ex- 
perimental runs. 

The approximate ranges over which the independent and dependent variables 
ranged in this series of experiments are listed in table I. The dc power input to 
the plasma ranged over at least six orders of magnitude, from less than l/lO of 
a watt to more than 150 kilowatts. At present, arcing within the vacuum tank and 
flashover of the external insulators limits the potential which can be applied to 
the electrode rings, and the total current and power drawn by the plasma are 
limited by uncooled sheet metal and glass viewports in the vacuum tank. 


C ORRELATION PROC EDURE 

The dependent variables chosen to be correlated include the ion number den- 
sity and kinetic temperature; and the plasma current flowing to the electrode 
rings, which is a dependent variable characteristic of the modified Penning dis- 
charge ion heating scheme used in the NASA Lewis Bumpy Torus. The plasma 
current is also a measure of the total particle loss rate from the confined plasma. 
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In the absence of an adequate understanding of the physical processes in the plas- 
ma, it was assumed that these three dependent variables could be correlated by 
power law relations of the form 


r _ c «11«S1 «41 ^51 

r p" C l V A P 0 B max n e 

(1) 

^22 ^22 ^42 a 52 

T i = c 2 VVVVaxV 

(2) 

. . c “13 “23,.“33 “43 
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Each dependent variable was assumed to be a power law function of the three 
variables which can be independently varied in the laboratory, the anode voltage 
in kilovolts, the background neutral gas pressure in the vacuum tank measured 
in microtorr, and the maximum magtt* tic field on the magnetic axis of the plasma 
in tesla. Also included were the electron number density in particles per cubic 
centimeter and plasma current in amperes, which were considered possible in- 
dependent variables of a physical theory describing the dependent variables. 
Exceptions were equations (1) and (3), where die dependent variable was not in- 
cluded on the right-hand side of the equation, that is, is identically equal to 
zero in equation (1), and a 53 was identically equal to zero in equation (3), 

Having chosen the power laws in equations (1) to (3), the possibility must be 
recognized that one or more of the parameters included on the right-hand side of 
these equations may not be an appropriate correlation parameter. A parameter 
may not be appropriate for at least two reasons; The parameter may be irrele- 
vant to the physical processes which determine the quantity on the left-hand side 
of the equation; or the parameter may be inappropriate because two or more of 
the variables on the right-hand side of the equation may be implicit functions of 
one another, or related to each other by physical processes which do not allow 
them to be treated as though they were independent variables. In *he second 
case, such implicit dependences among the variables violate the fundamental as- 
sumption that the variables being correlated are statistically independent. 

To allow for the possibility that one or more of the variables might be ir- 
relevant in equations (1) to (3), one or more of the exponents on the right-hand 
side of these equations were set equal to zero according to the scheme described 
in table II. In set no. 1, all of the variables in equations (1) to (3) were permitted 
to vary. In sets no. 2 to 9, one or more of the "independent" variables was omit- 
ted from the correlation, 
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The power law expressions in equations (1) to (3) were converted to simple 
linear polynomials by talcing the logarithms of both sides and treating the expo- 
nents as unknowns which were to be determined by obtaining the best fit 
consistent with the experimental data for the variables which they multiply. The 
best-fitting power law exponents rv^ were obtained by a statistical procedure 
which is described in reference 17. This procedure not only provided the best- 
fitting values of consistent with the experimental data of table III, but it 
also provided two quantities which are useful in specifying quantitatively how 
good is the fit to the experimental data. The first of these quantities is the RMS 
normalized residual, which is given by the square root of equation 4. 22 in ref- 
erence 17. For brevity, this RMS normalized residual will be called the "er- 
ror, * * and gives a measure of how far the experimental data points lie from the 
best-fitting pow ' law given in equations (1) to (3). The larger the error, the 
further away are the experimental points from the value predicted by the power 
law relationship. 

An example of the physical significance of the "error" may be seen in fig- 
ure 1, in which the measured value of the plasma current is plotted on the x-axis 
and the right-hand side of equation (1) for the best- fitting values of a^. is plot- 
ted on the y-axis for run series UY-VJ. If the correlation lav' were exact, all 
of the data would lie along the straight line with slope of 45° and the error would 
be zero. In this instance, the data for the low pressure mode, designated by the 
letter "A", has an error of 0. 221; and these data points lie very close to the 
line of agreement. However, the high pressure mode data, designated by the 
letter "B", has an error of 1.01, five times as high; and the spread of the data 
from the line of agreement is much greater. Naturally, one would like this error 
to be as small as possible. The smallness of the "error" was used as the pri- 
mary criterion for choosing among the nine sets of correlations listed in table II. 

The correlation procedure also produced a second quantity which is equal to 
the square root of the left-hand side of equation 4, 23 in reference 17, which we 
will call the " pseudo- variance, " This pseuCo- variance may be thought of as the 
error bars on the parameters • .. The pseudo- variance becomes smaller as the 
fit to the experimental data becomes better, and it becomes smaller as the popu- 
lation of data points which is used to determine the exponents a. j becomes 
larger. Large values of the pseudo- variance imply a small number of data points 
and/or a poor fit to the experimental data. 

EXPERIMENTAL DATA POPULATION 

The experimental data, which was correlated according to the scheme de- 
scribed above, is listed in table III. These data were taken over a period of 
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2 years. The first three run series used spectroscopic techniques to measure 
the electron temperature. The relative electron number density was obtained 
from the electron temperature by a procedure which is described in references 
3 and 5. Unfortunately, only relative and not absolute number densities could be 
measured by this technique. Run series YM-YQ, ZL-ZM, and ZNl-ZNfi used the 
RF emission diagnostic described by Gerdln in references 4 and 12 to measure 
the absolute plasma number density. 

The data shown on table III represents both helium and deuterium gas, and 
both positive and negative polarities of the midplane electrode rings. There is 
one run series, YC-YJ, in which a single electrode ring waB used instead of the 
12 used in all other run series. Unfortunately, the electron number density was 
not measured for this run series. As a result, it was not useful in deriving cor- 
relations for comparison with the case of 12 electrode rings. Such comparisons 
will be the subject of future research. 

DETERMINATION OF BEST FITTING SCALING LAWS 

The experimental data listed in table III were analyzed by a computer pro- 
gram which performed the statistical procedure described in reference 17 for 
one run series at a time. For each run series, separate correlations were ob- 
tained for the high and low pressure modes of operation. The run series listed 
in table III were not combined into one large data population for two reasons. 

The first was thav the computer memory capacity was inadequate to handle all 
of the data listed on table III. Even if this limitation could have been overcame, 
it was anticipated that the exponents of the scaling law or the coefficients Cj 
multiplying the right-hand side of equations (1) to (3), might vary with the mode 
of operation, the gas used, the polarity of the midplane electrode rings, or the 
number of electrode rings used. 

The first order of business in correlating the experimental data was to de- 
termine which of the nine combinations of parameters listed in table n best cor- 
related the available experimental data. All of the run series on table III were 
analyzed by the computer program, and correlations were obtained for equations 

(1) to (3) and for all nine sets of variables. The best- fitting correlation among 
the nine sets on table II had the following characteristics. 

(1) The error was minimized for the best-fitting combination of parameters. 

(2) The pseudo- variance of the exponents a:., was minimized for the best- 
fitting set of variables. 

(3) Five of the run series listed in table III were taken under identical condi- 
tions. The values of the best-fitting exponents should be the same for these data. 
For sets of parameters that were not the best- fitting ones, the values of 




could and did differ substantially anu i hese otherwise identical run series. 

For the best-fitting exponents, however, these five run series yielded nearly the 
same values of the exponents and the multiplying constants Cj. 

Inspection of this mass of information revealed that the sets of parameters 
with the wavy underline in table II were those which produced the best- fitting o'f- 
relations to the experimental data, That is, set 7 best fit tire current correia^on 
of equation (1), set 3 best fit the ion temperature correlation in equation (2), and 
set 6 best fit the electron number density correlation in equation (3), The best- 
fitting correlation laws therefore had the functional form listed in equations (4) 
to (6). 
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Some specific examples of the best-fitting correlation laws are shown in fig- 
ure 1 to 3 for one of the run series on table III. In these figures, the letter 
"A M denotes data taken in the low pressure mode of operation, and the letter 
"B" data from the high pressure mode. On the ordinates are plotted the right- 
hand side of equations (4) to (6), with the best-fitting values of the constant Cj 
and of the exponents for that particular run series. The experimentally mea- 
sured value of the correlated parameter is plotted on the x-axis. 

The characteristic data on figure 1 show that the data are reasonably well 
correlated over about three orders of magnitude in the measured plasma current, 
and there do not seem to be any systematic departures from the correlation laws 
either at the low or the high end of the operating range. Figure 2 shows the best- 
fitting correlations of the ion kinetic temperature for the same run series. The 
data scarcely range over an order of magnitude. This run series is early data in 
which the charge- exchange neutral detector was not at its maximum sensitivity, 
and in which the ion kinetic temperature was limited by the anode voltage avail- 
able from the power supply which was used at that time. Data from another run 
series, YT-YZ, are illustrated in figure 2(b) in which the ion kinetic temperature 
varies over an order of magnitude. The data for both modes of operation are well 
correlated and do not vary from the correlation law at either the low or the high 
end of the operating conditions. On figure 3 are shown the best-fitting correlations 
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of the electron number density for the same run series illustrated in figures 1 
and 2. In this case the correlations range over about two orders of magnitude 
and represent reasonably good fits to the data, although there appears to be dif- 
ficulty at the highest densities in figure 3 where the high pressure mode varies 
substantially from the best-fitting correlation. 

Since we are interested in extrapolating the results of the current experi- 
ment, it is important to determine if there are departures from the scaling laws 
n«ar the upper limit of the parameters being correlated. This was done by de- 
termining whether there were knees oi’ asymptotes in the plotted data on figures 
1 to 3, and for similar such graphs for the other run series listed in table III, 

In addition to looking for departures from the line of agreement, the data wore 
further analyzed by plotting the ratio of the right-hand side of equations (1) to 
(3) to the experimentally measured value of that parameter, as a function of the 
independent variables being correlated. This procedure is illustrated for fig- 
ure 3, in which the correlated number density has a suspiciously wide spread 
at the high end of its range, suggesting that the cor 'elation might break down at 
high densities. To investigate this behavior, figure 4 shows the ratio of the cor- 
related to the measured value of the density as a function of the three independent 
variables being correlated. In figure 4(a), there is no tendency, in either mode 
of operation, for the correlated number density to vary systematically with the 
neutral gas pressure, even though the data are scattered rather widely at the 
higher values of neutral gas pressure. In figure 4(b) this ratio is plotted as a 
function of plasma current, and again there is no systematic departure from a 
good correlation either at the low or high end of the plasma current. Figure 4(c) 
shows the ratio plotted as a function of anode voltage, and there may be a slight 
failure of the scaling law at high anode voltage. Such a departure from the 
scaling law, however, was not observed in the other run series which were cor- 
related. 


DETERMINATION OF SCALING LAW EXPONENTS 

After determining that the combinations of parameters underlined in table II 
produced the best- fitting correlations for the experimental data, the next objec- 
tive is to determine whether the best-fitting values of the exponents a'.j differ 
systematically with the conditions under which the data were taken, that is, 
whether the exponents differ with the type of gas used, the polarity of the elec- 
trode ring, or the mode of operation. The best-fitting exponents for each run 
series which could be correlated are shown in figures 5 to 7 along with the 
pseudo- variance, which provides error bars for the exponent. There did not 
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appear to bo any systematic dependence of '',o exponents on the type of gas used, 
the polarity of the electrode rings, or the mode of operation. There did not seem 
to be any systematic effect on the ion kinetic temperature correlation of the single 
case on flgui'e 7 in which only one electrode ring was used. 

The vertical lines in figure 5 to 7 represent an effort by the author to define 
a compromise exponent which would apply to the entire data population. These 
overall best-fitting value of are listed in table IV for the scaling laws given 
in equations (4) to (6). The total number of experimental runs represented by the 
data on figures 5 to 7 are summarized in the third column of table IV, and the 
constants which multiply equations (4) to (6) are listed under table IV. While 
there did not appear to be any systematic dependence of the exponents on the op- 
erating conditions listed, the constants Cj which multiply the scaling relation- 
ships did depend on theso external factors. 

The value of the constant C^, which determines the current draw by the 
plasma, depends both on the mode of operation and on the polarity of the elec- 
trode rings. 

The constant Cg, which multiplies equation (6) for the ion kinetic tempera- 
ture correlation, does not depend on the polarity of the electrode rings, the mode 
of operation, or the gas used. This is illustrated by data from reference 4, 
shown in figure 8, in which the ion kinetic temperatures were measured under 
identical conditions of gas pressure, magnetic field, and anode voltage with both 
positive and negative polarities applied to the midplane electrode rings. In fig- 
ures 8(a) and (b) the ion kinetic temperature and the ion heating efficiency do not 
depend on the electrode polarity. That the constant Cg is apparently the same 
for both deuterium and helium gas is somewhat surprising, although the theoret- 
ical expression for the ion kinetic temperature derived in reference 15 predicts 
that the ion kinet- c temperature is proportional only to the fourth root of the par- 
ticle mass ratio. The 16 percent difference that results from use of helium in- 
stead of deuterium would probably not be resolvable with the current data. 

The coefficient C g multiplying the density correlation differs only with the 
polarity of the midplane electrode rings, and is about a factor of 4 lower for nega- 
tive polarities than with positive electrode rings. It is interesting to note on fig- 
ure 7 that the exponents of the density correlation are the same for those run 
series in which the relative number density was measured by spectroscopic 
means, and for those in which the absolute number density was measured by the 
RF emission technique of Gerdin. This agreement of the scaling relationships 
for the two sets of density data reinforces Gerdin’ s conclusion (refs. 4 and 12) 
that the RF emission peak near the lower, hydrid frequency enabled him to mea- 
sure a quantity proportional to the plasma number density. On figure 9 are shown 
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some data from reference 4 in which Gerdin plotted the relative density from 
spectroscopic measurements on the ordinal and the absolute number density 
from the RF measurements on the abscissa. They tended to fall around a 45° 
line of agreement. Gerdin' s data were charactoidstic of an emitting region in 
the plasma, which may be a region of maximum density. The spectroscopic data, 
however, represented a time and space average across the plasma diameter, 

PHYSICAL IMPLICATIONS OF SCALING LAWS 

The best-fitting scaling laws may be written down by substituting the expo- 
nents from table IV in equations (4) to (6) to yield: 
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where the constants C^ to C 3 are given on table IV. When scaling the plasma 
current, it is often desirable to have only independent variables on the right- 
hand side. This may be obtained by substituting equation (9) into equation (7) to 
obtain 


amps 

where the constant is given by 
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4- 
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These scaling laws have some interesting implications. To the extent that 
these exponents represent a good compromise among the individual run series 
shown in figures 5 to 7, then this implies that the physical processes which gov- 
ern the four parameters correlated in equations (7) to (9) are the same regard- 
less of the mode of operation, the type of gas, or the polarity of the midplane 
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electrode rings. This similarity of the physical processes is particularly sur- * 
prising for the two modes of operation, In references 3 and 4 evidence was pre- 
sented which seemed to indicate that different physical processes may have been 
affecting the radial transport of particles in the two modes of operation. The 
above correlation laws appear to suggest that the physical processes occurring 
in the plasma are the same in the two modes, and the different currents repre- 
sented by the differing values of in the two modes of operation are merely 
the result of a geometric or boundary condition effect, rather than representing 
an entirely different physical process. 

The electron number density in equation (9) does not depend on the mode of ' , 
operation. The value of Cg is the same In both modes for the RF derived num- 
ber density measurements, but the relative number densities obtained spectro- 
scopically differ between one and two orders of magnitude for the two modes. 

The value of Cg in table IV is taken from the RF data. The reason why a sub- 
stantial modal density difference occurred in the spectroscopic data but not in 
the RF data must await determination by future measurements with a microwave 
interferometer, 

The relationship between the plasma current and the other variables shown 
in equation (7) may be understood on the following basis: The continuity equation 
for charged particles in the Bumpy Torus plasma may be written in the form 


dn 
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where the first term on the right-hand side is the volume production rate of 
charged particles in the plasma, and the second term is the volume rate of loss, 
where the plasma current ^ is the total rate of loss of particles from the 
plasma, and Vp is the plasma volume. In the steady state, one can set the 
left-hand side of equation (11) equal to zero and obtain 


1 p = Q V p n e n 0 <av) ne (12) 

Equation (12) demonstrates that when the plasma current is due to volume ioniza- 
tion, it is proportional to the product of the neutral and charged particle number 
density. This is very close to the empirical scaling law of equation (7). It is 
therefore possible to understand equation (7) as a consequence of the plasma 
current being created by volume ionization in the plasma, and not by conditions 
on the boundary, collection of parasitic currents, or by emission from one or 
another of the electrodes used. 
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Equation (8) for the ion kinetic temperature is not tiro functional dependence 
which one would expect from simplistic physical arguments. It has been shown 
in references 3, 4, 14, IB, and 10 that the ions are heated by ExB rotation in 
the crossed olectric and magnetic fields of tho plasma. On this basis, one might 
expect the ion kinetic tomporature to be given in terms of the E/B drift velocity 
as follows, 
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where the average olectric field is assumed to be the electrode potential divided 
by the plasma radius. The oversimplified model of equation (13) therefore pre- 
dicts an ion kinetic temperature proportional to the ion mass, proportional to the 
squai 3 of tire electrode potential, and inversely proportional to the square of the 
magt**tie field. 

!'n reference 15 a theory taking into account finite gyroradius effects, and 
which assumes that the ion mobility is higher than the electron mobility, yielded 
the theoretical Expression for the ion kinetic temperatime given in equation (14). 



eV 


(14) 


This theory predicts that the ion kinetic temperature is proportional to the 
fourth root of tire ion mass, to the fourtlr root of the electron number density, 
inversely proportional to the square root of the magnetic field, and directly pro- 
portional to anode voltage. 

The empirical scaling law in equation (8) is much closer to the theoretical 
expression derived for the modified Penning discharge in equation (14) than it is 
for the simple ExB drift model. The weak dependence on electron number den- 
sity predicted by equation (14) may be beyond the resolution of the current data 
and does not necessarily represent a discrepancy. The differences between the 
exponents on the electrode potential and magnetic field are substantial, and prob- 
ably represent a discrepancy larger than the uncertainty of the current data set. 


DISCUSSION 

Several parameters which are known to affect the plasma characteristics 
have not yet been optimized. These parameters include the number of electrode 
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rings used to generate the plasma; the alignment of the oloctrode rings with re- 
spect to the magnetic field drift surfaces; and the presence, direction, and mag- 
nitude of a weak vertical magnetic field imposed parallel to the major axis of the 
torus, Systematic variation of thoso threo parameters has been awaiting installa- 
tion and operation of a mlorowave interferometer system so that reliable absolute 
electron number density measurements could be made. If changing these threo 
parameters degraded the plasma characteristics and made the scaling laws dis- 
cussed above less favorable, one could leave them at thoir existing values with- 
out penalty. 

There is reason to suppose that the variation of one or more of these param- 
eters will substantially improve tho electron number density obtained. On fig- 
ure 10 is shown some spectroscopic data taken by R. W. Richardson which shows 
the relative electron number density as a function of the number of anode rings. 
The electron number density increases substantially as tho number of anode rings 
is decreased- tho low pressure mode of operation, the relative electron num- 
ber density iv eraasod by about a factor of 5 as the number of anode rings was re- 
duced from 12 to 6. In the high pressure mode of operation, the relative electron 
number density increased by about a factor of 10 as the number of anode rings 
was reduced from 12 to 3. If these trends are confirmed by subsequent data 
taken with a microwave interferometer, it may be possible to increase the con- 
stant which multiplies the scaling law of equation (10) for the electron number 
density by a factor of 5 to 10 if the optimum number of electrode rings are used. 

Figure 11 suggests what may be an effect on the ion kinetic temperature and 
ion heating efficiency of the number of anode rings used to generate the plasma. 
Shown on figure 11(a) is the ion kinetic temperature with 12 anode rings on the 
ordinate, and the ion kinetic temperature measured under the same conditions 
with only a single anode ring used to energize the plasma on the abscissa. There 
does not appear to be any significant difference in the ion kinetic temperatures in 
the low pressure mode of operation, but in the high pressure mode of operation 
12 electrode rings appear to heat the plasma more effectively than a single one. 
This same general trend for the same data population is exhibited in figure 11(b) 
which shows the ion heating efficiency for 12 and 1 anode rings. 

Thus, the available data seem to indicate that it may be possible to improve 
on the results implied by the scaling laws derived above, if the remaining experi- 
mental parameters are optimized. 
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SUMMARY 


Parametric variation of independent variables which may affect the charac- 
teristics of the NASA Lewis Bumpy Torus have identified those which have a sig- 
nificant effect on the plasma current, ion kinetic temperature and number density, 
and those which do not. Empirical scaling laws which relate these parameters to 
the midplane electrode voltage, magnetic field, and neutral gas pressure were 
obtained from a series of experimental runs containing several hundred data 
points and an order of magnitude or more variation in the latter variables. The 
best-fitting scaling laws were found to be as follows. 
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In addition to the continuous independent variables correlated above, additional 
parameters which affect the plasma properties were also varied including the 
polarity of the midplane electrode rings (which determines the direction of the 
rad'al electric field), the number of midplane electrode rings, the mode of 
plasma operation, and the method of measuring the plasma number density. 

The scaling laws derived from the data remained valid over an order of 
magnitude or more in the dependent parameters being correlated and showed no 
tendency to break down at the highest plasma •. „ ents and ion energies and den- 
sities which were attainable in the current ser,° i of experiments. The condi- 
tions at the upper limits were sufficient to exhibit small but measurable steady- 
state neutron production from DD reactions in a deuterium plasma. 
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TABLE II. - COMBINATIONS OF VARIABLES CORRELATED 


Set 

number 

Ip, function 1 

Tj, function 2 

i\ e , function 3 




“31 



“l2 

"22 

“32 

“42 

“52 

“13 

“23 

“33 

a 

CO 

“53 

1 

V 

— 

V 

V 


V 

V 

V 

V 

V 

V 

V 

V 

V 


2 

V 

— 

V 

V 

0 

V 

V 

0 

V 

V 

V 

0 

V 

V 


3 

V 

— 

V 

0 

0 

V 

JO 

- 

m v m 

V 

V 

V 

0 

V 

— 

4 

V 

— 

V 

0 

V 

V 

T* 

*v 1 

V 

’ v’ 

V 

0 

0 

V 

— 

5 

V 

— 

0 

V 

V 

V 

0 

V 

V 

0 

V 

V 

V 

0 

— 

6 

V 

— 

0 

A 

l / 

V 

V 

V 

V 

V 

0 

V 

JU, 


0 


7 





JL 

V 

0 

V 

0 

V 

0 

V 

V 

V 

— 

8 

0 

— 

V 

V 

V 

V 

V 

0 

0 

V 

0 

V 

0 

V 

- 

9 

0 

— 

V 

0 

V 

V ! 

o 

0 

0 

V 

0 

0 

V 

V 

— 


V variable to be correlated 
0 exponent set equal to zero 

— variable not in correlation 






PARAMETERS VARIED DURING EACH EXPERIMENTAL RUN SERIES 





























I 


I 


18 


TABLE IV. - PRELIMINARY SCALING LAWS FOR BUMPY TORUS 


Variable 

correlated 

Best 

fitting 

set 

number 

Total 

experi- 
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runs 

Constant 

Parameter exponents 
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Figure 1. - The best- fitting correlation of Eq. *4) for the 
plasma current, as a function of the measured plasma 
current tor run series UY-VJ. Letter "A" Is the low 
pressure mode of operation, and "B" is the high pres- 
sure mode of operation. 



(a) RUN SERIES UY-VJ. 

Figure ?. - Best fitting correlation of ion kinetic temperature 
from Eq. (3) as a function of the measured ton kinetic tern 
perature. The letter ''A" is the low pressure mode of oper 
atton and "B" designates the high pressure mode of oper- 
ation. 
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NIGH PRESSURE MODE 


MEASURED ION KINETIC TEMP. T1, eV 
lb) RUN SERIES YT-YZ. 
Eiqure 2. * Concluded. 
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Figure \ - The best-fitting correlation of Eq. (61 is 
i function of measured absolute electron number 
density, for run series UY-VJ. The letter "A" 
designates the k ># pressure mode of operation, 
and "8" the high pressure mode. 
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Figure 4. - Ratio of the best fitting correlation of Eq. i6i for 
the electron number density In run series UY-VJ, divided 
by the measured *lue of electron number density as a 
function of the parameters varied during this run series. 
The lettt. "A" designates the low pressure mode of oper 
ation, and "B" the high pressure mode of operation. 
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Figure 4. Continued. 
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RINGS POSITIVE, T, 4 , ?f tV 
tailON KINETIC TEMPERATURE. 

Figure 8. - Date from ref. 4 Indicating that 
the Ion kinetic temper etu re end ion heat- 
ing efficiency do not depend on electrode 
ring polarity or mode of operation. Data 
were taken for the same magnetic fields, 
electrode voltage and background neu 
tral gas pressures. 



ION HEATING EFF WITH ALL 12 EUCTRODE 
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IbllON HEATING EFFICIENT. . 

Figure 8. - Concluded. 
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Figure 9. Relationship between the relative 
ni'Ti oer density measured spectroscopically, 
on the y-axis, and the Asolute electron 
number density, measured by RF emission 
methods, on the x axis. The stralqht line 
with slope of 45 degrees Indicate that both 
diagnostic methods were measuring a qur. 
tity proportional to electron number density. 



Figure 10. * Ejtperi mental data taken from ref. 3 show- 
ing the relative electron number deputy obtained 
spectroscopically as a function of the number of anode 
rings for *he high and low pressure mode of operation. 
Bmax * 2. 4 tesla, and helium gas. 
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ill ION KINETIC TEMPERATURE WITH ALL 1? ELEC- 
mODE RINGS POSITIVE PLOTTED AGAINST THE 
ION KINETIC TEMPERATURE WITH ONLY ONE 
EUCTRODE RING ENERGIZED AND FOR THE SAME 
CONDITIONS OF MAGNETIC FIELD. BACKGROUND 
PRESSURE. AND ANOOE VOLTAGE. 

Figure 14. * The effect of the number of anode 
rings on Ion kinetic temperature and Ion heat- 
ing efficiency, taken from ref. 4. 
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ib) THE K)N HEATING EFFICIENCY WITH ALL 1? Et£C- 
mOOE RINGS POSITIVE COMPARED WITH THE ION 
HEATING EFFICIENCY WITH ONLY ONE ANODE RING 
ENERGIZED FOR THE SAME CONDITIONS OF BACK- 
GROUND PRESSURE. MAGNETIC FIELD. ANO AN- 
ODE VOLTAGE. 

Figure 11. - Concluded. 
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